1 99 6 BU - CCS - 950601 BROWN - HET - 999 hep - lat 9509012 MIT - CTP - 2461 Chronological Inversion Method for the Dirac Matrix in Hybrid Monte Carlo
نویسندگان
چکیده
In Hybrid Monte Carlo simulations for full QCD, the gauge fields evolve smoothly as a function of Molecular Dynamics time. Here we investigate improved methods of estimating the trial or starting solutions for the Dirac matrix inversion as superpositions of a chronological sequence of solutions in the recent past. By taking as the trial solution the vector which minimizes the residual in the linear space spanned by the past solutions, the number of conjugate gradient iterations per unit MD time is decreased by at least a factor of 2. Extensions of this basic approach to precondition the conjugate gradient iterations are also discussed.
منابع مشابه
Simulation of full QCD using overlap fermions
Most lattice Quantum Chromodynamics (QCD) formulations explicitly break chiral symmetry. This leads to unwanted mixing between operators of different chiral sectors, exceptional configurations when the simulated quark masses approach their physical values, and no clear definition of a topological index. The overlap Dirac operator solves these problems by satisfying a lattice variant of chiral s...
متن کاملNumerical methods for the QCD overlap operator IV: Hybrid Monte Carlo
The extreme computational costs of calculating the sign of the Wilson matrix within the overlap operator have so far prevented four dimensional dynamical overlap simulations on realistic lattice sizes, because the computational power required to invert the overlap operator, the time consuming part of the Hybrid Monte Carlo algorithm, is too high. In this series of papers we introduced the optim...
متن کاملThe overlap lattice Dirac operator and dynamical fermions
I show how to avoid a two level nested conjugate gradient procedure in the context of Hybrid Monte Carlo with the overlap fermionic action. The resulting procedure is quite similar to Hybrid Monte Carlo with domain wall fermions, but is more flexible and therefore has some potential worth exploring.
متن کاملM ar 1 99 6 Exact Local Bosonic Algorithm for Full QCD
We present an exact local bosonic algorithm for the simulation of dynamical fermions in lattice QCD. It is based on a non-hermitian polynomial approximation of the inverse of the quark matrix and a global Metropolis accept/reject correction of the systematic errors. We show that this algorithm is a real alternative to the Hybrid Monte Carlo algorithm.
متن کاملImplementing Hybrid Monte Carlo with stout-smeared chirally improved Dirac operators
We discuss our implementation of dynamical Ginsparg-Wilson type fermions using a stoutsmeared chirally improved Dirac operator. Such operators have been studied extensively in quenched calculations within the Bern-Graz-Regensburg (BGR) collaboration. Here we discuss the development and testing of the Hybrid Monte Carlo algorithm with this Dirac operator. We study the chiral properties of this o...
متن کامل